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Abstract 

Our previous results [J.Phys.: Condens. Matter 14 (2002) 13777] dealing with the 
analytical solution of the two-dimensional (2-D) Anderson localization problem due 
to disorder is generalized for anisotropic systems (two different hopping matrix 
elements in transverse directions) . We discuss the mathematical nature of the metal- 
insulator phase transition which occurs in the 2-D case, in contrast to the 1-D case, 
where such a phase transition does not occur. In anisotropic systems two localization 
lengths arise instead of one length only. 
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1 Introduction 



Anderson localization [T] remains one of the main problems in the physics 
of disordered systems (see e.g. the review articles [2|3|4] ). In the series of our 
previous papers [5]l6f7] we presented an exact analytic solution to this problem. 
By the exact solution we mean the calculation of the phase diagram for the 
metal-insulator system. We have been able to solve the two dimensional (2-D) 
problem [5]. We have shown then that the phase of delocalized states exists 
for a non-interacting electron system. The main aim of the paper [6] was 
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the generalization of the results to the case of higher dimensional spaces (N- 
D). In paper [7j we discussed the mathematical details of the new analytical 
approach for calculating the phase-diagram. An exact solution is only possible 
for the conventional Anderson model: the tight-binding approximation with 
diagonal disorder, where on-site potentials are independently and identically 
distributed. 

It is well known that exact results in the field of phase transitions (the metal- 
insulator transition is a particular case of a phase transition) are exceedingly 
rare [gpj . This is why any extension of the applicability range of analytical 
methods in this field is of great interest. In this paper we extend our ap- 
proach [5|6|7p0] for anisotropic media [TT|12lll3|[Ti] , but else remaining in the 
framework of a conventional Anderson model. 

Incorporation of an anisotropy into the tight-binding approximation with di- 
agonal disorder is also methodologically valuable. Before the exact solution 
was obtained for D > 1 [5]|6|7] , analytical methods concentrated on the D = 1 
case p!5|IT6] . However, the specific topology of 1-D systems does not permit 
extension of the results to higher dimensions. This is true in particular for the 
Ising model, where an exact analytical solution for the 1-D case has nothing to 
do with the 2-D solution obtained by Onsager [8lf9] . As it is well known for the 
Anderson localization problem, all states in the 1-D system are localized (i.e. 
there is no metal-insulator transition). This is a particular result of a general 
theory [2j that no phase transitions are possible in 1-D systems with short- 
range interactions. In fact, phase transitions (e.g. in the Ising model [511^ or 
the Anderson localization problem j5]) are observed only starting with D = 2. 

It should be realized that the approximate methods are also of a limited use 
here. In particular, traditional perturbation theory works perfectly in the 1-D 
case [T^, which is well demonstrated by the analysis of the exact solution in 
Ref . piTT] . Random potentials can be treated in the 1-D Anderson localization 
problem as a small parameter, which is used in series expansions of physical 
quantities. However, this approach fails [9] for systems with phase transitions 
in D > 1, since physical quantities here are no longer described by analytical 
functions and series expansions. The same is valid for the Lyapunov exponent 
7 (which is the inverse of the localization length ^ in the Anderson problem 

0). 

Of particular interest is the understanding of the mathematical nature of the 
phase transition: how does the analytical character of the exact solution for 
the 1-D problem changes to a non-analytical character of the exact solution 
for the 2-D case? For the discrete spatial dimensions (1-D or 2-D) a simple 
comparison of the two relevant solutions does not help us. However, a ^^contin- 
uous" treatment of the spatial dimension for anisotropic systems, performed 
in this paper, provides much more insight. We have put the word continuous 
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in quotation marks because we have in the present case of anisotropic systems 
an interesting possibihty of a transition from a 2-D system to a 1-D system, 
by considering the hmit k, ^ for parameter of anisotropy k. 

The paper is organized as follows. In Section 2 we discuss the basic equations 
for the isotropic problem [5], which are generalized there for the anisotropic 
case. Section 3 presents the main results for the anisotropic problem. We 
demonstrate how the study of the hmiting case of a strong anisotropy permits 
us to establish a relation between 1-D and 2-D cases. 



2 Recursion relation and the filter function 

2.1 Isotropic system 

Let us start with the Schrodinger equation for the isotropic system (the lattice 
constant and the hopping matrix element are set equal to unity) 

n,m- V / 

The on-site potentials en,m are independently and identically distributed with 
existing first two moments, {en,m) = and (sl^ .^ = cr^. 

Eq.([T]) can be written as the recursion 

+ C.1pn,m, (2) 

where the operator C acts on the index m and is defined by the relation 

(3) 

m'=±l 



The standard initial (boundary) conditions are ipo^m = and ipi^m = C(m, 
respectively. 

The existence of some fundamental 1-D numerical series /i„ (n = 0, 1, . . . , oo), 
the so-called system function or filter, was proved rigorously in Ref. [5)16117] . A 
study of the asymptotic behaviour of this series allows to define uniquely the 
phase diagram of the system. Namely, the series /i„ is bounded 

\hn\ < OO, (4) 
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for the delocalized states, but increases without bound 

hm \hn\ oo, (5) 



n— >oo 



for the locahzed states. The actual numerical values are not important in 
this respect. The determination of the existence region of the bounded series, 
eq.(|l]), permits to find also that for the delocalized states, i.e. to obtain the 
system's phase diagram. In its turn, the asymptotics of the series (jS]), /i„ oc 
exp(27n), for the localized states allows us to extract the Lyapunov exponent 
7, or the localization length ,^ = 1/7. 

It is more convenient to determine the phase diagram and the Lyapunov ex- 
ponent 7 by using the complex variable z (Z-transform) 

= £ ^- (6) 

n=0 ^ 



In this case the physical problem of the localized/delocalized states is reduced 
to the mathematical problem of the search for the poles of the complex variable 
function H{z) [51I61I7] . which is also called the filter. Eq.(l5]) means in this case 
that H{z) has the pole at z = exp(27). 

For the isotropic system the filter can be found exactly pi6|7] : 



rr 1 / N . Z + 1 / rf/c 



w 



2 i^+^r 



where 

C{k) = E ~2cos{k) (9) 



corresponds to the operator C and arises due to use of the Fourier transform. 

Note that solution of the Schrodinger equation ([2]) depends on chosen initial 
conditions, Vo.m = and ipi^rn = ctm- In contrast, as follows from the defini- 
tion, eq.([7]), the filter H{z) is independent of the initial conditions (i.e. the 
field variables am), and thus is not the functional of the wave function ipn,m- 
Following [5|6|7] . let us briefly explain the mathematical nature of the filter 
H{z). According to the definition, eq.Q, the operator C acts on the wave 
functions. The expression eq. ([7]) is a functional of C and a disorder parameter 
a. The analysis [5115117] shows that the function H{z) is also an operator which 
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acts, however, not on wave functions but certain mathematical objects called 
signals. Only these signals are functionals of the wave functions. 

The relation between the filter H(z) and signals (input signal S^^\z) and 
output signal S{z)), 



S{z) = H{z)S^^\z), (10) 
(0) 1 z+1 J \a{k)\^dk 

— 7r 

has a certain similarity with the transition to the operator formalism in quan- 
tum mechanics. Note that the boundary conditions (field or Fourier trans- 
form a{k)) influence only functions S^^\z) which are independent of the pa- 
rameter a. In other words, functions S^^\z) correspond to solutions in a com- 
pletely ordered system. Introduction of disorder, cr > 0, transforms the initial 
solution S'^'^\z) into S'(z), where the operator H{z) describes this transforma- 
tion. 

The concept of the filter function is a general and abstract description of the 
problem of localization. Instead of analyzing wave functions, i.e. different sig- 
nals, it is sufficient to analyse properties of the fundamental localization opera- 
tor H{z) by means of the theory for functions of complex variables. The signal 
concept (input and output signals) and the filter function are key elements in 
the theory of signals used earlier |5|i6|17j . In these papers the connection be- 
tween the Anderson localization problem and signal theory is discussed and 
the most important concept of the proposed method, the filter H{z), is defined. 
The knowledge of papers [5)16117] presenting the main formalism is prerequisite 
for understanding the present paper. 

Let us discuss equation ([5]). The localized states correspond to values of 7 > 
0, i.e. a divergence of the filter function /i„ (or signals). This divergence is 
mathematically similar to the divergence of averages in the diffusion motion 
[7]. That is, transition from delocalized to localized states can be treated 
as a generalized diffusion with a noise-induced first-order phase transition. 
The generalized diffusion arises due to the instability of a fundamental mode 
corresponding to averages of wave functions (correlators and signals). 



2.2 Anisotropic system 



The Schrodinger equation for the 2-D anisotropic medium 

) + t2{lpn,m+l ) = {E - en,m)'^n,m (12) 
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contains two hopping matrix elements: ti,t2- Let us use the units, where 
max{ti,t2} = 1, and min{ti,t2} = >^ with < k < 1. The parameter k 
characterizes the system's anisotropy. For example, the bandwidth is \E\ < 

= 2(l + /€). 

Eq. f|T2l) can also be written as a recursion relation. Let us assume e'^^^ = 
£n,m/ti and define the operator C through the relation 

^ ^ m'=±l 



As a result, we obtain the recursion 



functionally similar to eq.(l2]). The equation for the filter function H{z) trans- 
forms, respectively, to 



where {a'Y = {{e'Y) = a' /ti 



E t2, 



C\k) = ^2cos(A;) (16) 

ti ti 



In an anisotropic system there are two distinct directions. We give them dif- 
ferent names in order to distinguish them later. The choice ti = 1, ^2 = 
corresponds to the direction of strong binding, while ti = ^2 = 1 to the 
direction of weak binding. Note that the former case (strong binding) should 
only be considered for establishing the correspondence with the exact 1-D so- 
lution [51IT7] in the limit k — >^ 0. Indeed, for k = one gets ti = 1, t2 = and 
the operator C = E. Eq. fll4l) transforms, respectively, into the ensemble of 
similar equations characterized by a neutral index m, 



each of which is equivalent to the 1-D equation 

1pn+l = -£ni^n " + Etpn- (18) 
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2.3 Useful relation 



Let us use the relation 

I } dk 1 



277 J C±2cos(A;) i^/A^'' 



(19) 



which holds provided the imaginary part of the complex variable C is non- 
negative, ImC, > (C is defined in the upper semi-plane). 

The integrand f|T5l) can be presented as 

1 11 1 



w2 - (/:'(A;))2 2w ^^x; - a{k) w + C'{k) 



and then eq. (fT9|) can be used, provided that the parameter w is also defined 
in the upper semi-plane, Imw > 0. 

The transformation from the complex variable z to w {conformal mapping) 
gives a parametric presentation of the filter function H{z) = T-C{w) [6] conve- 
nient for analysis. Since the relation between the variables z and w is nonlinear, 
this relation is not unique. 



2.4 Conform transformation and phase index P 



If the complex variable w = u + iv is defined in the upper semi-plane, w > 0, 
the relation 

w = P{z'/^ + (21) 



is valid [5], where the parameter P = ±1. As shown below, P serves as an 
index numbering different-type solutions: localized (P = +1) and delocalized 
(P = —1), i.e. it characterizes the phase state of a system. 

For P = +1 eq. (pT]) characterizes the conform transformation of the outer 
region of the unit circle, \z\ > 1, onto the upper w-semi-plane. Similarly, for 
P = — 1 the conform transformation of the inner part of the unit circle, \z\ < 1, 
occurs. 

The inverse transformation is described by the relations 
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z = — 



IH h P — V4 

2 2 



— w'^, 



(22) 



provided 



z+1 



P 



w 



z — 1 iy^A — w"^ 
2.5 Anderson localization and first-order phase transition 

As was shown [SE], for arbitrary dimension D the phase of focahzed states 
(the solution with P = +1) is found for an arbitrary disorder a. In contrast, 
the defocahzed phase (P = —1) arises only for D > 2. It exists within the 
limited energy band and below the critical disorder parameter a. 

A comparison of the 1-D and 2-D cases shows their qualitative differences. In 
1-D even infinitesimal disorder destroys the localized states; their existence 
range is degenerate into a single point cr = [5|ITF] . Similarly to many other 
problems of phase transitions [8], the (metal- insulator) phase transition in 1- 
D is impossible. In this respect the delocalized states at a = 0, which are 
unstable under any perturbation, cannot be treated as a special phase. It is 
more reasonable to treat these states as the limiting case of the localized states 
which only exist for 1-D. As a consequence, the physical characteristics of the 
single-phase system such as the Lyapunov exponents, 7 = j{a,E), are an- 
alytical functions of the energy E and the disorder parameter a. Moreover, 
the perturbation theory also applies here: 7 = 7(0", E) can be expanded in 
a series of a small parameter jSj. The 1-D system is characterized by a 
continuum transition due to the lack of disorder, a 0. The limiting value 
\im„^Q'y{a, E) = 7(0, i?) has obviously to coincide with that for the delocal- 
ized states, 7 = 0. In other words. 



The 2-D case differs qualitatively from the 1-D case. First of all, in the existence 
range of the second solution (P = —1) the first solution (P = -|-1) is also 
always present. In other words, two phases co-exist 05] ■ Therefore, the metal- 
insulator transition has to be interpreted as a. first-order phase transition. Two 
fundamental points should be stressed here, (i) In any system with dimension 
larger or equal to two and with a phase transition, the physical properties are 
no longer analytical functions of their variables. That is, such properties as the 
Lyapunov exponents 7(0", E) can manifest peculiarities, e.g. step-like changes. 



7(0,P) = 



(24) 



holds. 
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and thus the expansion in small parameter o"^ no longer holds here [5]. (ii) 
The first-order phase transitions reveal additional peculiarities. The coexisting 
phases (in our case, the delocalized states with 7 = and localized states with 
7(0", E) 7^ 0) have nothing in common, and thus there is no reason to expect 
that eq. fl24|) holds. That is, in a system with first-order phase transition 

7(0,E)^0, (25) 



which illustrates once more the above-discussed non-analytical behaviour of 
the physical properties. 

As was recently noted [5], 7(0", £") in the Anderson localization problem can 
be treated as the long-range order parameter. Indeed, the two different phases 
reveal different 7 parameters, 7 = and 7 7^ 0, respectively. The first-order 
transition is always characterized by the step-like change of the long-range 
parameter [9J. In this respect, eq. (l25|) being the particular case of a more 
general relation 7 7^ for the phase of delocalized states, seems quite self- 
evident, despite possible exceptions discussed below. 

Therefore, the analysis of the limiting transition to 1-D in the anisotropic 
system permits, on the one hand, to understand the mechanism for the change 
from non-analytical to analytical solutions for the physical properties, and, on 
the other hand, gives an additional check of our results [5|l6] since the difference 
between eqs.f lMj) and fl25l) is nontrivial. 



3 Results 



3.1 General remarks 



Using the above-obtained relations, one gets the following parametric presen- 
tation of the filter function 

where 

Q,{w,E) = ^-^w + u^. (27) 

Let us consider now the properties of the filter function H{w) = Ti^iw) at 
P = —1. As was shown in Ref. [516], the filter analysis permits to obtain the 
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phase diagram, i.e. to determine the existence range of the delocahzed states. 
For this one has to determine the poles of H{z), i.e. to find the roots of the 
equation H~^{z) = 0. In the parametric presentation for P = — 1 we seek 
roots of 

nz\w) = 0, (28) 



or 



E ; . = 1- (29) 



The pecuharity of the case P = — 1 [5|6] is that its physical interpretation is 
possible either in the absence of the poles (which is the case for D > 2 [6J), or 
if the poles lie on the unit circle \z\ = 1 [marginal stability), i.e. 

z = exp{±i2<^) (30) 



(pairs of complex values), ip G (0,7r). The letter case applies to 2-D isotropic 
systems [5] . It can also be shown using the method in ref . [6j that taking into 
account the anisotropy does not affect eq. (!30|) (we omit here the proof). As a 
result, in the general case of an anisotropic system it follows from eq. fl2ip that 

w = —2 cos{(p), (31) 



i.e. the w parameter is real. This obviously simplifies the solution of eq. (l29!) . 
The equation 

nz\w = 0) = (32) 



determines the threshold disorder a = (Jq{E) which defines the limits of exis- 
tence of the delocalized phase [5|6]. 

It is easy to conclude that the real roots w of eq. fl28p satisfy the condition 
\Qy{w,E)\ < 2, which is equivalent to 

ti\w\ + \E\<2t2. (33) 



Simultaneously < \w\ < 2 has to hold. 
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Consider now the filter function T-C-^{w) for P = +1. This function has always 
the only pole [5] 



z = exp{2j), (34) 

which lies on the real axis with 7 > 0, where 7 is interpreted as the Lyapunov 
exponent. In the parametric presentation the parameter w is also real, 

w = 2 cosh(7). (35) 

To find the Lyapunov exponent, one has to determine the root of the equation 
n~\w) = 0, (36) 



or 

,2 



+ ^=^=— — ] = 1- (37) 



2tit2Vw^ - 4 jQl(w,E)-A jQl{w,E)-A 



The parameter w sought for satisfies the natural condition of positive argu- 
ments of all roots of eg. (1371). so the simple condition 

tiw - \E\ > 2t2 (38) 
has to be fulfilled, along with w > 2. 



3.2 Isotropic system 



Let us summarize for further analysis the main results for an isotropic system. 



3.2.1 The case of P = —1 (delocalized states). 

In an isotropic system the real roots w have to fulfill < < 2 [5J. This is 
possible only in the limited energy range l^'l < Eq = 2, otherwise one obtains 
complex roots with an unclear interpretation. That is, the delocalized states 
exist only within a certain energy range around the zone center. In this energy 
range an increase of the disorder parameter a shifts the poles in the parametric 
presentation from \ w\ = Wq = 2 — \E\ < 2 towards the limiting value of w = 0. 
For further disorder increases, the real roots disappear and thus the formal 
solutions of eq. (j28il permit no physical interpretation. 
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Using eq. (l32l) at ti = ^2 = 1, one gets [5] 



ao{E)=2{l-E'lElf'\ 



(39) 



3.2.2 The case P = +1 (localized states). 

Eg. (1551) determines the existence of real roots w > Uq = 2 + \E\ oi eq. flTn) 
OS]- The eq. flS^ holds here with 7 = 'y{a,E). In the parametric presentation 
the value of w = uq corresponds to the absence of disorder, cr = 0. In this case 



This result confirms eq. fl25|) . The continuity, eq. (1241) . is found only at the band 
center, E = 0, when 7(0, 0) = 0. In this case one obtains for the Lyapunov 
exponent 



E.g. for cr — s> one has from eq. (HT]) 7(0", 0) oc a. It follows from this that 
the function 7(0", E) cannot be represented as a series in powers of cr^, i.e. 
perturbation theory is not applicable to the 2-D Anderson problem. 

3.3 Strong binding direction 

Let us now assume that ti = 1, t2 = k,- As was mentioned in section fI72[ this 
case in the limit of k ^ corresponds to the transition to the 1-D case. 

3.3.1 The case P = — 1 (delocalized states). 
In this case eq. (l33l) transforms into 



Along with the relation < \w\ < 2 one obtains again that the real roots 
exist in the energy range of \E\ < Ei = 2k,. A detailed study shows that this 
situation is similar to that for the isotropic system, with a marginal stability 
of the filter function. 

The threshold disorder value reads now 




(40) 



7(cr,0) 



sinh \-). 



(41) 



w\ + \E\ < 2k. 



(42) 
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ai{E)=2^{l-E''/ElY'\ 



(43) 



Therefore, on can see a monotonic decrease of both energy half- width Ei (x k 
(where delocahzed states exist) and the threshold disorder cri(£') oc y/H in 
the limit k ^ in the strong binding direction. For k = the phase of 
delocalized states loses its existence range, which corresponds to the well- 
known 1-D result. 

3.3.2 The case of P = +1 (localized states). 

Eg. (1551) transforms into w > 2k + \E\. Together with w >2 one gets 

w>Uo = max{2,2K+\E\}. (44) 



For the lack of disorder, a — > 0, w ^ Uq. 

It is easy to find that at \E\ < E2 = 2{1 — k) the parameter mq = 2. Respec- 
tively, in this energy range for a = the Lyapunov exponent 7(0, E) = 0. At 
small a the Lyapunov exponent is proportional to cr^: 7(0", E) ~ /3cr^. Analysis 
of eq. flS^ yields 

^{E, k) = ][ I ^ + , ^ ] (45) 

4 J{E^a. + E) {E2 + E) J{E„,a. - E) {E2 - E) 



In the 1-D limit, k — 0, one gets E^ax = 2, £'2 = 2, \E\ < 2, respectively 
/?(^'0) = ^^. (46) 



As one can expect, this corresponds to the exact 1-D result (in the limit 
of small a): 

-f{a,E) ^ a^iA- E"^). (47) 



In the energy range E2 < \E\ < Emax Uq = 2k + \E\ > 2. Respectively, the 
value 

\E\ 

7(0,E) = cosh"i(/t + ^) (48) 



is nonzero. 

Therefore, in the isotropic system with /t = 1 due to the presence of the phase 
transition eq. lHOl) confirms completely the general conclusion, eq. (l25il . The 
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only exception is the band center E = 0, where 7(0, 0) = holds, and there 
is a correspondence with the main relation for the 1-D system - eq.( l24p . As 
anisotropy arises, k > 0, this point moves into the energy range \E\ < E2, 
where eq.f l24p holds. Moreover, within this range 7(0", E) can be expanded into 
a series in a small parameter a^, similarly to the 1-D case. Outside this energy 
range eq. fl25l) remains valid. As we have demonstrated, at k = our equations 
transform into those obtained earlier for the 1-D system [5|T7] . 

3.4 Weak binding direction 

Let us assume now ti = k,, t2 = I. 

3.4-1 The case P = — 1 (delocalized states). 
In this case the following relations 

k\w\ + \E\<2 (49) 

and ()<\w\<2 hold. 

Real roots of equation 7iZ^{w) = exist in the parameter range < |w| < 
min{2, (2 — \E\)/k} for energy \E\ < Eq = 2. The threshold disorder magni- 
tude, where the delocalized states still exist, remains to be defined by eq.f lS^ : 

a2iE) = 2^(1 - E''/ElY'\ (50) 

3.4.2 The case P = +1 (localized states) 
Eq. (!38l) corresponds to 

w>uq = {2+\E\)/k, (51) 

whereas the second condition w > 2 is fulfilled automatically. As a ^ 0, 
w uq> 2 holds. Therefore eq. fl25l) holds in the whole energy range. 

3.5 Phase diagram 

The delocalized state in the anisotropic system can be discussed only provided 
delocalization occurs simultaneously in two dimension. It is easy to see that 
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decisive point here are relations obtained for the strong binding direction. 
Therefore, the phase of delocahzed states exists in the range \E\ < Ei = 2k, 
provided for a given energy the disorder does not exceed the critical value, 
eq. (H5]l . A decrease of the anisotropy parameter k leads to a monotonic de- 
crease of both the energy range and the threshold disorder. 



4 Conclusions 

The phase of localized states in the anisotropic system can be characterized by 
two Lyapunov exponents 7(0", E) and, respectively, two localization lengths. 
The Lyapunov exponent in the strong binding direction is always smaller then 
that in the weak binding direction. In other words, the localization length in 
the strong binding direction is always larger then that in the perpendicular 
direction. These two localization lengths reveal different asymptotic behaviour 
as (T ^ 0. In the limit of strong anisotropy, k —>■ 0, only the localization length 
in the strong binding direction serves as an analogue to the single localization 
length in the 1-D system. 
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